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^ . Abstract 

O 

o 

Q"^ I In one-loop string effective action, we study a generality of non-singular cos- 

■*-«» . niological solutions found in the isotropic and homogeneous case. We discuss 

O^' Bianchi I and IX type spacetimes. We find that nonsingular solutions still 

zTp.' exist in Bianchi I model around nonsingular flat Friedmann solutions. On the 

other hand, we cannot find any nonsingular solutions in Bianchi IX model. 
^ , Non-existence of nonsingular Bianchi IX universe may be consistent with the 

j^ ■ analysis by Kawai, Sakagami and Soda, i.e. the tensor mode perturbations 

against nonsingular flat Friedmann universe are unstable, because Bianchi IX 
model can be regarded as a closed Friedmann universe with a single grav- 
itational wave. So based on these facts, we may conclude the nonsingular 
universe is found in isotropic case is generally unstable, a singularity avoid- 
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ance may not work in the present model. 
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I. INTRODUCTION 

The initial singularity is one of the most serious problems in the Big Bang universe. 
Even the inflationary universe model, which resolves many difficulties in the early universe, 
cannot provide us an avoidance of initial singularity. Quantum effect might resolve it as a 
scenario in quantum cosmology, but quantum gravity is not yet completed. To construct a 
theory of quantum gravity, a unification of fundamental interactions may be one of the most 
promising ways. Among many attempts of such a unification, a superstring model would 
be the best candidate []I|. In fact, the origin of a black hole entropy may be understood by 
a string theory. Hence, the initial singularity problem might also be solved in the context 
of a string theory, which may provide us to manage the physics at Planck scale. However a 
full theory have not yet developed. We may so far not be able to discuss the early universe 
in a string theory as it is, except for some restricted models [0. Therefore many efforts 
to find nonsingular universe are mainly based on a low energy effective field theory of a 
superstring model. Such an effective field theory will break down beyond the Planck energy 
scale. However, if a superstring theory is truly the theory of everything, its effective theory 
may reveal some important aspects about gravity. In particular, if it has a property of 
singularity avoidance, it would be a nice evidence for a superstring. 

Based on an effective theory, there is one interesting approach called pre-big-bang uni- 
verse model 0. This model is based on T-duality, which gives some relation between large 
and small scales. Then, assuming that a scale factor of the Universe has such a duality [Q , 
we find a cosmological solution which consists of two distinct and disconnected branches. 
One of the branches(t > 0) corresponds to the expanding Friedmann universe, and other 
branch(t < 0) gives another expanding universe, ending with a singularity at t = 0. In pre- 
big-bang scenario, however, those two disconnected branches are assumed to be connected 
without a singularity on account of some unknown stringy effect, which is not included in 
the lowest effective action. 

There are several works to study whether two branches can be connected without a 



singularity. At the tree-level in a superstring action, however, a possibility of classical 
branch-changing solutions is excluded 0. Only quantum effect may be a hope to connect 
two branches. In fact, taking into one- loop contributions to a superstring effective action 
with dilaton and modulus fields |0, Antoniadis, Rizos and Tamvakis found an existence of 
nonsingular solution with spatially flat background, with the same action, then, Easther and 
one of the present authors also showed a nonsingular closed universe 0. Although those 
solutions do not help the pre-big-bang scenario, they may be still interesting because of 
singularity avoidance. 

A question about such nonsingular solutions, however, is whether those are generic or not. 
Although we find nonsingular solutions for some finite range of initial data, the spacetime 
is assumed to be isotropic and homogeneous 0. The universe, however, may begin with 
an anisotropic and/or inhomogeneous geometry. Therefore, we have to study whether the 
present nonsingular solution is generic or not, or is stable or unstable. Recently, Kawai, 
Sakagami and Soda analyzed stability of flat nonsingular solution against perturbations, and 
found that a tensor mode is unstable. However, since a closed nonsingular universe bounces 
in a finite time, it is not clear whether it is unstable against tensor mode perturbations. 
In order to clarify such a problem and to study generality of nonsingular solutions, in this 
paper, we will analyze two types of Bianchi models; Bianchi I and IX model. 

This paper is organized as follows. In the next section, we introduce the basic equations 
for Bianchi I and IX models. The numerical results are presented in Sec. Ill, and conclusions 
and discussions follow. We adopt the metric signature (—,+,+,+) and units of c = SvrG = 1. 

II. BASIC EQUATIONS 



We take the following one-loop effective action p] - pA\ 
S = J d'xV^l^R - \^^^^'^ - ^V^aVV - ^-H^'-'H,,, + ^[Ae* - (5e(o-)]i?^5 

+ (Higher Curvature Terms) , (2.1) 
where R, (j) and a are the scalar curvature, the dilaton, and the modulus field, respectively. 



^GB = R — ^R^'^R^iv + R^ RfiuaP (2.2) 

is the Gauss-Bonnet term, and H is the antisymmetric tensor field. The coefficient A is 
positive definite and determined by the inverse string tension a' and string couphng g. 
The coefficient S depends on the relative numbers of chiral, vector and spin-| massless 
supermultiplets and is proportional to the four- dimensional trace anomaly of the N = 2 
sector. The 6 would be either positive or negative. The function ^(cr) is defined as 

a^) = ln[2eV(^e'^)] (2.3) 



by the Dedekind t] function [12 



ry(r) = gTi n(l-g2^), g = e^"^ (2.4) 

n=l 

The first derivative of ^ with respect to a is 

•^ „=i ■■- e 

which can be approximated by sinh a as 

ea(^)~-Ysinhcr. (2.6) 

As was shown in |^, Eq. ( |2.6| ) is very close to the exact function ( |2.5| ). Then we will use 
Eg. (^.61) in our analysis just for simplicity, although we have also checked our results by use 
of the exact function. We also introduce the following function /(0, a) for convenience: 

fi<f>.^) = ^[e''-mcr)], (2.7) 

where 6 = 6/X. We set H = and ignore higher curvature terms than second order. 

With the same model, cosmological solutions have been studied from a view point of 
the initial singularity problem. Antoniadis, Rizos and Tamvakis analyzed a spatially fiat 
Friedmann model and found nonsingular solution 0. Easther and Maeda extended their 
analysis to a closed Friedmann model and also showed nonsingular solution ^. Since we wish 
to examine generality of those nonsingular solutions, we extend their works to anisotropic 



spacetimes. Here we study only Bianchi I and IX models because those spacetimes include 
a flat and a closed Friedmann models. 

It maybe convenient to rescale time and spatial coordinates by A as t = t/vA, x* = 
xYvA. Hereafter we drop a bar for convenience. Taking a variation of the action, we obtain 
the basic equations. With those basic equations, for Bianchi I and IX models, we can assume 
the following diagonal metric form: 



ds' = -de + e-'''e'^^^u'u:^, 



(2i 



(3+ + v^/5_ 

A.J = (3+- 73/3. 

-2(3^ 

In both cases, we also introduce scale factors a{t),b{t), and c(t) as 



gP{t) ^ g-n+/3++v^/3_ ^q(t) ^ g-n+/3+-x/3/3_ ^r(t) ^ g-f^-2/34 



(2.9) 



The basic equations obtained are divided into two groups: 
(1) The dynamical equations for the metric p, q and r, the dilaton field 0, and the modulus 
field a : 



(1 + 8f/)(g + f) + (1 + 8g/)(f + r^) + (i + 8/) 



qr + -f/i 



-'■12 , 3 2 



(1 + 8f/)(p + p2) + (1 + 8pf){r + r^) + (i + 8/) [rp + U2] + ^0' + ^a' = 0, 
(1 + 8qf){p + p') + (1 + 8pf) {q + q') + (1 + 8/) [pq + f/3] + ^0' + ^^^ = 0, 



+ (p + g + r) 



-i?: 



.dfp2 
25/02 



GB; 



a + (p + g + r)cj = ^g^^GB, 



(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 



where an overdots denotes a differentiation with respect to t, and the Gauss-Bonnet term 
Rgb is given as 



Rgb - 8 



111 

[ji + p^] [qr + -Ui] + [q + q^] [rp + -f/2] + [r + r^] [m + -f/3] 



+-p -7: \- -^q -7\ \- -r -7: 7;P(l[-7: h 



2 dp 2 dq 2 dr 2 dp dq dr 
2 dq dr dp 2 dr dp dq J 



(2) The constraint equation : 



1 • 3 

pq + qr + rp + 2Apqrf + ApfUi + 4g/f/2 + ArfU^ — -0^ o"^ 

+ ^(f/i + t/2 + f/3) = 0. (2.16) 

Here the functions Ui,Ui, Ui are defined by 
(A) Bianchi I model : 



Ui = U2 = U3 = (2.17) 

(B) Bianchi IX model : 

f/2 = e-2" + e-^P - e"^" + ^(e^^'^-P-^) + e^^^-"-'''^ - Se^^^-'-^)) 

f/3 = e-2p + e-^*? - e-2^' + l(e2(p-'?-'^) + g^C?— p) _ Se'^'^-^-^)). (2.18) 

The basic equations above are the five second-order differential equations with one con- 
straint equation. For p{t) = q{t) = r{t), i.e. the case of isotropic and homogeneous model, 
the equations are reduced to the cases of a flat universe studied by Antoniadis, Rizos and 
Tamvakis , and of a close universe by Easther and Maeda [0] . 

III. NUMERICAL RESULTS 

We have examined the case of 5 < because in the isotropic and homogeneous case, 
nonsingular cosmological solutions are found only for 6 < |^, |0|. We solve the basic 
equations numerically. 



Since five second-order derivatives p,q,r,(f),a in the basic equations (|2.10|) - (|2.14|) are 



coupled, we have to make an inverse transformation as follows: Defining a vector x 
{p, q, r, 0, a), the basic equations ( p.lO|) -( p.l4| ) are written in the matrix form as 



Zx = y, (3.1) 

where 5x5 matrix Z = Z{j>, q, r, 0, a, p, q, r, 0, &) and vector y = y{p, q, r, 0, a, p, q, r, 0, a) 
are known explicitly from the basic equations. 
Unless A = detZ vanishes, we have 

X = Z-'y, (3.2) 



then we can solve the basic equations (|3.2| ), giving initial data oi p,q,r,(l),a,p,q,r,(l),a. 
However, when A vanishes, then we cannot proceed our numerical calculations further. Such 
a end point, which may appear in the evolution of the Universe, seems to be a spacetime 
singularity, however more detail analysis will be required as we will show later. 

Although we have used scale factors p, q and r for the basic equations, we shall describe 
our results by two anisotropy variables /?+ and /?_. 

A. Bianchi I Type Case 

First we will show the results in Bianchi I model. We choose S = —AS/tc. We introduce 
scale factors as 

a{t) = e^W = e-^+'3++^'3-, 6(t) = e''^*) = 6"^+^+-^'^-, c(t) = e'^^*) = e-''-^f+ . (3.3) 

Without loss of generality, we can set ^o = P± = 0, i.e. ao = bo = cq = 1. Here the 
subscript denotes an initial value of the variables. We define the initial time t = Iq = 0. 
We have to give initial data of ^o, i3±q, 0o, 00, ctq and do, which must satisfy one constraint 



equation ( |2.16| ). Then we have five independent initial values. Since we are interested in 
whether nonsingular solutions found in the isotropic case are generic or not, we shall set up 
the initial data in the anisotropic case around the isotropic ones. 
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In the isotropic and homogeneous case, we find nonsingular solutions for some finite 
parameter range of initial data. In Fig. 1, fixing CIq = — 0.1,cro = we show the range of 
initial data of 0o and (pQ, which gives for nonsingular solutions (shown by a circle), do is 



determined by the constraint equation (|2.16|). In Fig. 1, if 0o is efficiently small, then in 



the case of 0o < there always exist nonsingular solutions. We think because if < and 
101 ^ 1, then e*^ ~ so singularity avoidance is almost independent of 0. 

To give initial data Qq, 0o, 0o and ctq with an anisotropy, taking one nonsingular isotropic 
solution, including anisotropy, i.e. /3+q,/?-_q(7^ 0), and solving for do by the constraint 
equation we have set up our initial data. 

We find nonsingular solutions in Bianchi I model near the isotropic nonsingular solution. 
We show one example in Fig. 2, where we have set CIq = —0.1, 0o = 0o = 0, ctq = and 
P+Q = 0.05, /9_Q = 0.0, do = 0.17320508075689, which has been determined by the constraint 
equation (do = 0.200 for the isotropic case). We also show a singular solution in Fig. 3. The 
volume e^ vanishes at t ^ —4.98 and t ^ 8.35, resulting in a Big Bang type singularity. 

Changing the anisotropic parameters /5+q, /?_q, we search for nonsingular solution in order 
to find how large anisotropy is possible to allow nonsingular solutions. The result is shown 
in Fig. 4 for Qq = —0.1, 0o = 0o = 0, a^ = 0. We find that if anisotropy is enough large 
at the initial stage, then the spacetime evolves into a singularity. The boundary between 
non-singular and singular solutions in (/5+q, /5_|_Q)-plane is almost a circle. Ar/r < O(10~^). 
Whether spacetime will evolve into a singularity or avoid it, of course, strongly depends on 
initial parameters which we have set {CIq, 0o, 0o, o"o)- As we move these initial parameters to 
the boundary values in the isotropic case beyond which no non-singular solution is found, the 
radius of the boundary in (/3+Q,/?+Q)-plane decreases and if we set initial parameters fairly 
close to the boundary values, then the range for non-singular solution eventually disappears. 
Then in this case the singularity avoidance does no longer work. 

We can conclude that for the Bianchi I type anisotropy, a singularity avoidance is still 
generic. Because the range of non-singular solutions is finite and not small. 

However, we have the stability analysis done by Kawai, Sakagami and Soda Jl^ [|l^ ||16 



Their analysis is based on the same action (but only a dilaton field) and fiat Friedmann 
background case. They found there exists instability in the tensor mode. Then in the case 
of the fiat universe the universe is unstable. This result may not change even in Bianchi I 
background universe. 

B. Bianchi IX Type Case 

In the same way in the Bianchi I type case, we will analyze Bianchi type IX. We introduce 
scale factors as 

a{t) = 2ePW, b{t) = 2e«W, c = 2e^W. (3.4) 

The factor 2 is definitions of invariant basis and the scale factor in a closed Friedmann 
model. As the same as Bianchi I model, we first research for closed nonsingular Friedmann 
solutions. We search for non-singular Bianchi IX solutions. 

First setting up the initial anisotropic data around those in a closed nonsingular solution, 
we show an isotropic non-singular solution in Fig. 5. We choose a negative value of 6 as 
S = -48/7r, and set a(0) = 6(0) = c(0) = e^-^^^^, CIq = -0.01, 0o = -3, 0o = -5 x 10-^ 
o"o = which was found by Easter and Maeda. We can reproduce one of their solutions. 

An anisotropy is induced just as the same as Bianchi I model, we show our numerical 
results in Fig. 6 and Fig. 7. As for anisotropy we set /3+q = 0.01 and /3_q = 0.0 in Fig.5, 
and /9+Q = 10"^"' and /?_g = in Fig. 6, respectively. 

In the evolution of the universe, A(= detZ) eventually vanishes. We always find a 
singularity at a certain finite time, so our numerical analysis have to break down. This type 
of breakdown appears even if an initial anisotropy is quite small, that is, both /3+q and /9_q 
are less than 10~^^. As shown in Fig. 6 and Fig. 7,1 = R^^'^^ R^jjap suddenly diverges near 
the breakdown point. 

We show the time evolution of (3+ and (3+ in Fig. 8 and Fig. 9, where (3+q = 0.01, 
/9_Q = 0.0 and /5+q = 10^"*^^, /3_q = 0, respectively. /5+, /3_ or both will also increase very 
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rapidly near the breakdown point. Because of the exponential growth of / = R^'^°'^ R^uap 
and (3±, we believe it is a curvature singularity. Since the volume factor e^ does not vanish 
there, it is not a Big Bang type singularity. The Universe find a singularity with a finite 
volume. 

We have searched the wide range of parameters where Easther and Maeda found non- 
singular closed universe solutions, but we could not find any nonsingular solution even if 
the anisotropy is extremely small. In all cases we examined I = R^^'^^ R^j^uap and (3± show a 
sharp increase near the point where A vanishes. Then we would conclude there is no non- 
singular solution in anisotropic Bianchi IX type models and all solutions except for exact 
isotropic case eventually evolve into a curvature singularity. 

Our results may be consistent with the stability analysis by Kawai, Sakagami and Soda 



T^ |T5| |16| . King showed that Bianchi IX type model can be regarded as a closed Friedmann 



background with a single gravitational wave of a fixed wave number k = yQ/ S where S is 



the radius of three-sphere |T^. This means that Bianchi IX model is the same as the closed 
Universe with a nonlinear tensor perturbation. Then, if the tensor mode is unstable even 
for a closed universe we can understand our results. In fact since we know the wave number 
A;, we can estimate the timescale of the tensor mode instability from their analysis and 
compare it with our numerical results. We find that our breakdown time is the same order 
of instability timescale, i.e. the order of our breakdown timescale is almost (0.1 ~ l)i/|(5|, 
which depends on initial parameters, and the instability timescale is ~ Wl*^!- 

IV. CONCLUSION REMARKS AND DISCUSSIONS 

In this paper we have examined generality of previously found cosmological nonsingular 
solutions in the heterotic superstring effective action in orbifold compactifications with one- 
loop correction. 

In Bianchi I type case, many nonsingular solutions are found around nonsingular flat 
Friedmann solutions without fine-tuned initial conditions. On the other hand, in Bianchi 
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IX type case, we cannot find any nonsingular solution even if anisotropy is quite small. 
In this case, anisotropy grows exponentially, and our numerical analysis breaks down at a 
certain point where anisotropy will diverge. At this point, both j3±, which represents the 
anisotropy, and / = R'^'^'^^ Rfiua/B seems to diverge, then the universe may evolve into a 
curvature singularity. The volume factor remains a finite value. 

Our result is consistent with the stability analysis by Kawai, Sakagami and Soda, who 
found tensor mode instability in fiat Friedmann background. Because Bianchi IX type model 
can be regarded as a closed Friedmann background with a single gravitational wave of a fixed 
wave number. 

Since the nonsingular fiat Friedmann model is unstable against tensor perturbations, we 
may expect that nonsingular Bianchi I model is not generic as well. Therefore, we may 
conclude the nonsingular universes found in isotropic cases are not generic and a singularity 
avoidance may not work in the present model. 

We have argued only first order expansion term of inverse string tension a' and one-loop 
correction and excluded antisymmetric tensor (a axion field) which might exist in the early 
universe. For more rigid consideration we may need to study based on a full superstring 
theory, though such a full theory have not yet been developed into a point where we can 
deal with cosmology. Then we hope to find some effects of a superstring theory by which 
Rt^'^^l^ R^^^^ ^ill j^ot diverge in some wide range of initial values. It might be given by higher 
curvature terms in the effective action. 
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Figure Captions 

FIG. 1: We show the fate of the flat Friedmann Universe in terms of initial 

values, 00 and <po in the cases of the flat Universe, o means regular 

solutions, while x means that the Universe evolve into a singularity. 

For A, the solution seems to be singular, although we could not 

confirm it because we need more CPU time. 
FIG. 2: We show one nonsingular solution. We choose a negative value of 

6 as 6 = —48/ 71. We have set Qq = —0.1, 0o = 0o = 0, ao = 

and /3+0 = 0.05, (3-^ = 0.0. &o = 0.17320508075689 has been 

determined by the constraint equation {&q = 0.200 for the isotropic 

case). We show the scale factor a in (a), dilaton field in (6), 

modulus field a in (c) and / = R'^'^°'^ R^ua/s in (d), respectively. 

The solid line, dashed line and dotted line represent scale factors a, 

b and c in the anisotropic case, respectively. The dash-dotted line 

represents a scale factor in the isotropic case. 

FIG. 3: We show one singular solution. We choose a negative value of S 

as (5 = — 48/7r. We have set Qq = —0.1, 0o = 0o = 0, (Tq = 

and /9+Q = 0.1, /3_q = 0.0. &o = 0.00 has been determined by the 

constraint equation {&q = 0.200 for the isotropic case). We show 

the scale factor a in (a), dilaton field in (6), modulus field a in (c) 

and I = R^^^^ Rfj^yaii in {d), respectively. The solid line, dashed line 

and dotted line represent scale factors a, h and c in the anisotropic 

case, respectively. The dash-dotted line represents a scale factor in 

the isotropic case. 
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FIG. 4: Changing the anisotropic parameters (3+Qand(3-Q, we search nonsin- 

gular solutions. For 5 = — 48/7r, Qo = —0.1 and ctq = 0, we show the 
boundary on the Z^+q and /9_q plane, beyond which no non-singular 
solution is found. The boundary is almost a circle. The solid line, 
dashed line and dotted line represent the boundary in the cases of 

00 = 00 = 0, of 00 = 1-5,00 = 0.20 and of 0o = 3.0, 0o = 0.40, 
respectively. 

Fig. 5: We show a nonsingular solution in a closed Friedmann universe. 

We choose a negative value oi 6 as 6 = — 48/7r. We set ao = e^'^^^^, 
Qq = —0.01, 00 = —3, 00 = — 5 X 10^"^, (To = 0. We show the scale 
factor a in (a), dilaton field in (6), modulus field a in (c) and 
/ = R'^""^^ R^uap in {d)-, respectively. 

Fig. 6: We show a singular solution in Bianchi IX type model. We set 

the same initial values in Fig. 5, 5^ = — 48/7r, ao = 6o = Co = 
e3-25i4^ r^o = -0.01, 00 = -3, 00 = -5 X 10-^ ao = 0, except for 
anisotropy, which is chosen as (3+q = 0.01 and (3-^ = 0.0. We show 
the scale factors a, b and c in (a), dilaton field in (b), modulus 
field (7 in (c) and / = R'^'^°'^ Rfj,uai3 in (d), respectively. We find that 

1 grows almost exponentially. 

Fig. 7: We show a singular solution in Bianchi IX type model. We set the 

same initial values as those in Fig. 6, except for anisotropy i.e. /9+q = 
lO""*^^ and /9_Q = 0. We show scale factors a, b and c in (a) and (b), 
dilaton field in (c), modulus field a in (d) and / = R^'^°'^ R^ua/s in 
(e), respectively. We find that / grows almost exponentially. 

Fig. 8: We show the time evolution of /5+ and /5+. We set initial parameters 

as those in Fig. 6. The anisotropy /3_|_ also grows exponentially, as 
I diverges. 
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Fig. 9: We show the time evolution of (3+ and (3+ for the solution given in 

Fig. 7. The anisotropy j3+ also grows exponentially, as / diverges. 
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